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Linear transport through a single-walled carbon nanotube ring, pierced by a magnetic field and 
capacitively coupled to a gate voltage source, is investigated starting from a model of interacting 
p z -electrons. The dc-conductance, calculated in the limit of weak tunneling between the ring and 
the leads, displays a periodic resonance pattern determined by the interplay between Coulomb 
interactions and quantum interference phenomena. Coulomb blockade effects are manifested in the 
absence of resonances for any applied flux in some gate voltage regions; the periodicity as a function 
of the applied flux can be smaller or larger than a flux quantum depending on the nanotube band 
mismatch. 
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I. INTRODUCTION 

Mesoscopic rings threaded by a magnetic field rep- 
resent an important tool for the investigation of quan- 
tum interference phenomena. The archetype example is 
the well known Aharonov-Bohm effect 1 , where the con- 
ductance of a clean ring exhibits a periodicity of one 
flux quantum $o = h/e. However, impurities 2 - and 
interactions 3 - can change this periodicity. In particular 
when transport through a one-dimensional ring is con- 
sidered, interactions lead to spin-charge separation, such 
that the most important contribution to the conductance 
arises when both charge and spin excitations propagate 
from one contact to the other arriving at the drain at the 
same time 3 -. Additionally, the dc-conductance of an inter- 
acting one-dimensional ring shows Coulomb oscillations 
with peak positions depending on the applied magnetic 
field and interaction strength^ 5 -. 

Among quasi-one dimensional systems, single-walled 
carbon nanotubes 6 (SWNTs) have been proved to be 
extremely interesting to probe electron-electron correla- 
tion effects. Luttinger liquid behavior, leading to power- 
law dependence of various quantities, has been predicted 
theoretically- 7 -* 8 - and observed experimentall y 9 ' 10 : 11 ' 12 in 
long, straight nanotubes. 

Moreover, as typical of low-dimensional systems, short, 
straight carbon nanotubes weakly attached to leads ex- 
hibit Coulomb blockade at low temperatures? . In metal- 
lic SWNTs two bands cross at the Fermi energy. To- 
gether with the spin degree this leads to the formation of 
electron shells, each accommodating up to four electrons. 
As a result, a characteristic even-odd^ or fourfol d 14 ' 15 ' 16 
periodicity of the Coulomb diamond size as a function 
of the gate voltage is found. Recently, spin-orbit effects 
in carbon nanotube quantum dots have been observed as 
welli 7 -. While the Coulomb blockade can be explained 
merely by the ground state properties of a SWNT, the 
determination of the current at higher bias voltages re- 
quires the inclusion of transitions of the system to elec- 




FIG. 1: Top view of a SWNT ring contacted to extended left 
and right leads and threaded by a magnetic flux 0. The ring 
is also capacitively coupled to a gate voltage V g plane situated 
beneath the ring (not shown). 



tronic excitations. Ini 8 - a mean-field treatment of the 
electron-electron interactions has been invoked to calcu- 
late the energy spectrum, while in 19 a bosonization ap- 
proach, valid for nanotubes with moderate-to-large diam- 
eters (0 > 1.5nm), has been used. Within the bosoniza- 
tion approach the fermionic ground state as well as the 
fermionic and bosonic excitations can be calculated. For 
small diameter nanotubes, short range interactions lead 
to pronounced exchange effect s 18 ' 20 which result in an 
experimentally detectable 16 singlet-triplet splitting. 

Generically, nanotubes have linear or curved shape. In- 
dividual circular single- wall carbon nanotubes have been 
observed i n 21 ' 22 ' 23 . However, these systems have been 
poorly experimentally investigated so far. Also from the 
theoretical point of view, only few works address the 
properties of toroidal nanotube s 24 ' 25 ' 26 ' 27 ' 28 ' 29 ' 30 . More- 
over, except for the studie s 25 ' 26 , on the persistent cur- 
rent and on the conductance of toroidal SWNTs, re- 
spectively, all the remaining theoretical works neglect 



electron-electron interactions effects. However, as for the 
case of straight SWNTs discussed above, electron correla- 
tion effects are expected to crucially influence the energy 
spectrum and transport properties of toroidal SWNTs. 
Indeed, ir>2£ it is found that for interacting SWNTs rings 
the persistent current pattern corresponds to the con- 
stant interaction model, with a fine structure stemming 
from exchange correlations. Results for the conductance 
of a SWNT ring weakly contacted to rings have been 
presented so far only in the short report, where a con- 
ductance resonance pattern depending on the interaction 
strength is reported. Moreover, a detailed derivation of 
the conductance formula for SWNT rings is missing in 
the short report^. In this work we generalize the analysis 
of^ on the conductance of interacting spinless electrons in 
a one-dimensional ring to the case of a three-dimensional 
toroidal, metallic SWNT at low energies. To this ex- 
tent we start from a model of interacting p z electrons on 
a graphene lattice at low energies and impose periodic 
boundary conditions along the nanotube circumference 
and twisting boundary conditions along the tube axis. At 
low energies only the lowest transverse energy sub-bands 
of the ring contribute to transport, so that the problem 
becomes effectively one-dimensional in momentum space. 
The three-dimensional shape of the SWNT orbitals in 
real space, however, crucially determines the final con- 
ductance formula. Indeed, in contrast to2&, we find the 
absence of interference terms between anti-clockwise and 
clockwise circulating electrons in the conductance for- 
mula. This is due to the localized character of the p z 
orbitals and the fact that for a realistic nanotube the 
contacts are extended. Additionally, we predict an eight 
electron periodicity of the conductance resonance pat- 
tern as a function of the applied gate voltage. Coulomb 
blockade effects are clearly visible in that for some gate 
voltage ranges a resonance condition is not met for any 
value of the applied flux. The resonance pattern is also 
periodic as a function of the applied magnetic field, with 
a periodicity which can be larger or smaller than one flux 
quantum for nanotubes with a band mismatch. 



II. THE TOTAL HAMILTONIAN 

We consider a ring made of a toroidal metallic single- 
wall carbon nanotube coupled to a source and drain elec- 
trode via tunneling contacts, cf. figure [TJ The torus 
is capacitively connected to a gate electrode with gate 
voltage V g , which changes the chemical potential of the 
ring. Further, a magnetic flux <fi threads the center of the 
torus. We wish to study how the two parameters V g and 
(j) influence the linear conductance of the system. The 
model Hamiltonian reads 



H = H r 



ing 



H L + H R + H T + H e 



(1) 



where -ffring, which also includes the parameters V g and 
<fi, describes the physics of the isolated SWNT ring and 
will be discussed in the next section. The second and 
third term refer to the metallic left and right contacts, de- 
scribed here as Fermi gases of non-interacting electrons. 
They read {a = L/R) H a = E^e^c^c^, where 
e a (g) is the energy dispersion relation of the lead a and 
c a qa is an operator annihilating an electron with wave 
vector q and spin a. 

The term Ht is the Hamiltonian describing tunneling be- 
tween the ring and the leads. Therefore it consists of two 
terms, Ht — Htl + Htr, and reads 



Hn 



EE 

x—R,L a 



d 3 r(T a (r)¥ a (r)$ aa (r) + h.c), (2) 



where ^(r) and $ (TQ (r) = Y^ff&tffflc&a are electron op- 
erators in the dot and in the leads, respectively, and T a {f) 
describes the generally position dependent transparency 
of the tunneling contact at lead a. In the following we 
choose the representation where f — (x, r±), with x being 
directed along the tube axis, while r± is the coordinate 
on the nanotube cross-section. With L being the circum- 
ference of the SWNT ring, the contacts are positioned in 
the region about x — and x = L/2, see Fig. [TJ Finally, 
H ext accounts for the energy dependence of the system 
on the external voltage sources controlling the chemical 
potential in the leads: H ext = -eJ2 a V a J2q„ c aqa c aqcr- 



The manuscript is organized as follows. The total 
Hamiltonian and its low energy spectrum are discussed 
in Sees. QTJ and IIII1 Specifically, the low energy Hamil- 
tonian of the non-interacting system possesses two linear 
branches, corresponding to clockwise and anti-clockwise 
motion, crossing at the two non-equivalent Fermi points 
of the graphene lattice. By inclusion of the dominant 
forward scattering processes only, the interacting Hamil- 
tonian is diagonalized exactly by standard bosonization 
techniques^. In IIVI the conductance formula is derived, 
while Sec. [V]is dedicated to the evaluation of the Green's 
functions of the interacting SWNT ring. Finally the con- 
ductance resonance pattern is discussed in Sec. IVI1 where 
also conclusions are drawn. 



III. LOW ENERGY DESCRIPTION OF 
METALLIC SWNT RINGS 

In this section we derive the low energy Hamiltonian 



#ring(0) = #1 



kin 



v.. 



ff ga te of the SWNT ring and 



diagonalize it. The ring Hamiltonian includes a kinetic 
term, the electron-electron interactions as well as the ef- 
fects of a capacitively applied gate voltage V g and of a 
magnetic flux </>. We assume a metallic SWNT. Hence, 
depending on the diameter of the SWNT, "low ener- 
gies" means an energy range of the order of I eV around 
the Fermi energy, where the dispersion relation for the 
non-interacting electrons near the Fermi points is lin- 
ear and only the two lowest subbands touching at the 



Fermi points can be considered. The effects of the mag- 
netic field are included by introducing twisted boundary 
conditions (TBC) 32 . Curvature and Zeeman effects are 
neglected here. Indeed the Zeeman splitting yields a con- 
tribution inversely proportional to the the square of the 
ring radius^ and is relevant only for very small rings. 
Finally, the linear dispersion relation around the Fermi 
points allows bosonization 31 of the interacting Hamilto- 
nian and its successive diagonalization when only the for- 
ward scattering part of the Coulomb interaction is in- 
cluded. The latter approximation is justified for SWNTs 
with medium-to-large cross-sections (0 > 1.5nm)2£. 



A. Twisted boundary conditions 

The band structure of non-interacting electrons in a 
SWNT ring is conveniently derived from the band struc- 
ture of the p z electrons in a graphene lattice. Since 
each unit cell of the graphene lattice contains two car- 
bon atoms, there are a valence and a conduction band 
touching at the corner points of the first Brillouin zone. 
Only two of these Fermi points, ±Kq = ±— 4=— e x , are 
independent. Since SWNTs are graphene sheets rolled 
up into a cylinder, we can obtain the electronic proper- 
ties of a SWNT by imposing quantization of the wave 
vector k = kj_ij_ + kueu around the tube waist, i.e., 
denoting the SWNT circumference with L±, one finds 
k± = j^-m, m — 0, ±1, ±2 . . . , which leads to the for- 
mation of several subbands labelled by m. We consider 
here only SWNTs of the armchair type, e x = eii, which 
are all metallic. In this case, only the two sub-bands 
touching at the Fermi points ±K~oe x are relevant. 

Figure[2]shows the linear valence and conduction bands 
of an armchair SWNT. Corresponding to the direction of 
motion of the electrons, there are two branches b = R, L 
with positive and negative slope, respectively. The linear 
dispersion relation reads 



s(b, k) = hvpsgn(b)n, 



(3) 



where k measures the distance between k\ \ and the Fermi 



points F = ±Kq, i.e., k 



II 



F, and vf is the Fermi 



velocity, vp = 8.1 x 10 5 m/s. Moreover, we use the con- 
vention b = r/l = ±1. The corresponding Bloch waves 
are of the form 



fFb K {r) 




(4) 



YjfpF^^W-R-fp). 



Here, Nl is the number of lattice points in the nanotube 
lattice. The index p = 1,2 refers to the two graphene 
sublattices and f p points from a lattice point R to a car- 
bon atom of sublattice p. For an armchair SWNT, the 
parameters f p pb are found to be f\Fb = 1/V2, /2F6 = 
— 1/\/2 sgn(Ffe). Finally, we have the p z electron orbital 

n. 




FIG. 2: Low energy spectrum of a non-interacting SWNT 
ring threaded by a magnetic flux <f>. The underlying graphene 
structure is reflected in the two Fermi points at ±A"o. Notice 
the misalignment </teo, with <f> — 4- + 5, of the energy levels 
between the left (L) and right (R) branches due to the intrin- 
sic mismatch 5 and the applied flux <j>. The flux quantum is 
denoted <j>o- 



As we consider a ring, we can impose in the absence of 
a magnetic field periodic boundary conditions along the 
direction parallel to the tube axis, yielding fc|| = ^m^, 
rnfc G Z. As the Fermi points ±Kq are not necessarily an 
integer multiple of this step size, we introduce a mismatch 
parameter S, S £ [—0.5,0.5), so that we can write K — 
=£-(n — 8), neZ. We still have to take account of the 
magnetic field. After travelling one time around the ring, 
an electron picks up a phase 27r^-, where <f>o — - is the 
flux quantum. This fact can be taken into account by 
introducing twisted boundary conditions (TBC), \&(r + 

Le x ) — e 1T ~*a^{f) , instead of periodic ones 3 ^,. Under 
TBC the allowed values of the momentum parallel to the 
tube axis k are 

2tt / <j) . 



2tt 

T 



(n K + 4>) 



n K € Z 

m K G Z. 



(5) 



Thus the electron operator ^ including the electron spin 
a =|, I reads 



a Fbon 

>(r)c F baK = \fL^2ipFb{r)i>Fba(x), 



e (fiFb( 



Fbcrfi 



Fba 



(G) 



where cpban annihilates an electron in state \Fban) and 
i/iFbaix) is a slowly varying ID operator. Later on the 



bosonisability of ipFba{x) will be of importance to evalu- 
ate the interference properties of the interacting SWNT 
ring. 



B. Kinetic energy and gate Hamiltonian 



The kinetic energy Hamiltonian in second quantization 
can in principle be read off from the dispersion relation 
for the low energy states (j3|): 

ffkin = fWF Yl S & n ( b ) KC Fba K C Fba K - ( 7 ) 

Fban 

Inserting the K-quantization ([5]), introducing the number 

counting operator Af F ba = E K c FbaK C Fba K and using tne 
abbrevation Eq = hvp^- we get 



flki. 



£ o E s § n ( & ) I E m " ( 

Fba \ K 



FbaK C FbuK + V^Fbcr 



(8) 



Gate voltage effects are included in the Hamiltonian 

■ffgato = -MgA/" = -eaV gl where TV = J^Fba^Fba is tne 
total electron operator of the SWNT, p g is the SWNT 
chemical potential and a a conversion factor. 



C. Interaction Hamiltonian 

A general form of the Hamiltonian V c - C describing in- 
teractions between the electrons is 



K. 



^ 



(9) 



x / d 3 r d 3 r , Vt(r^vL(r , )U(7 f -F , )V ( .,(i' , )i& a (f), 



where U{r — r ') is the possibly screened Coulomb poten- 
tial. By expressing the 3D electron operator in terms of 
the ID one, cf. Eq. ([6]), and integrating over the com- 
ponents of f and r' perpendicular to the tube axis one 
obtains an effective one-dimensional interaction. As dis- 
cussed ini, this yields in general interlattice and intralat- 
tice interactions. In2£ it is shown that when the tube 
diameter is large enough, the distinction between intra 
and inter-lattice interactions is no longer relevant, and 
only forward scattering processes, where the number of 
electrons in each branch remains constant, are relevant. 
Introducing the ID electron density operator 

PFba( x ) = ^Fba( x H 'Fba( x ) and retaining only for- 
ward scattering processes we obtain for the interaction 
Hamiltonian 



ft- - IEEE 

FF' bb' aa' 
X pFb<7(x)V (x,x')p F 'b'a'(x'), (10) 



L pL 

dx I dx' 

i) Jo 



with the effective ID potential Vq(x,x') defined as 
V (x,x') = ^J2 J d 2 r ± J d 2 r' ± x 



R.R' 



x \Ii{r-R)\ 2 U{r~r')\Il{r' - R 1 ) 



D. Bosonization 



(11) 



The problem of a one-dimensional system of interact- 
ing electrons is often conveniently formulated in terms 
of bosonic operators. A didactic overview can be found 
ir*2L. To this extent we assume a linear dispersion rela- 
tion over all values of the momentum k. This is justified 
by the long-range character of the interaction, leading to 
a momentum cut-off still lying in the low-energy range 
around the Fermi energy. The states with negative en- 
ergy are assumed to be filled in the ground state, the 
Fermi sea, of the system. As in^i, we now introduce new 
bosonic operators 



?Faq 



1 f E K C Fra K ' 

InJ 1 E« C Fla 



FraK+qi 



K^FltTK+qi 



q>0 

q<o, 



(12) 



with q = ^-nqjiiq G Z, and obeying the canonical 
bosonic commutation relations 



\ b F 1 a 1 q 1 ^ b F 2 a2q2^ ~ ° F 1 F 2 ""112 «9l92 > i b F 1 a 1 q 1 ' b F 2 a 2 q 2 J ~~ "' 



(13) 



In terms of the bosonic operators (TT2")) the kinetic energy 
Hamiltonian (TSl) assumes the form 



Hv 



£ Q^2^2\ n l\ h Faq h Faq 
Fa q^O 

|W1 



t() 



E 

Fba 



' Fb 



^ + 0sgn(6)AA J 



Fba 



(14) 



For the bosonization of the interaction part (jlQ|) , we de- 
compose the density operator pFba into its Fourier com- 
ponents, which are proportional to the bosonic operators 



P F ba( X ) = i>Fba( X )^Fba( X ) = -M 



Fba 



(15) 



E. I\n I f P is Sn(b)qx, , — jsgn(6)gxi t 

\J\""l\\y u Fasgn(b)q^ L U Fasgn 



(b)q ■ 



q>0 



By inserting (fTS")) into (JTUJ) we get 

V = \w m N 2 + E E E 5W9-9KI x 

Fa F'a' q=£0 
X ( b Faq + b Fa-q) ( b F>a>-q + b F>a>q) , 



with the total number of electrons Af = J^Fba-^Fba in 
the ring. The effective interaction Vq(x,x') is absorbed 



into 



L i-L 



dx / dx'e^ qx+q ' x 'W {x,x'). (16) 



W aa> 

/o Jo 

This is nonzero only for q' = —q and yields 

,L/2 

W q _ q =L- dycos(qy)V (\y\). (17) 

J-L/2 

The parameter Woo can be identified with the charging 
energy E c responsible for Coulomb blockade. To pro- 
ceed, it is convenient to introduce linear combinations of 
the bosonic operators (fT2|) associated to spin, charge and 
orbital degrees of freedom: 

b ®cq ■= -^{bK a ]q + b K „lq + b-K a ]q + b-K lq), 



?Gcq 



(bK„lq + b Kolq - b-K 1q ~ b-K lq), 



... • 7;( b K ]q - bK iq +b-K U ~ b-K a [q), 



■esq ■■= -j{b KaU -b Kaiq -b- K ^ q + b- Kaiq ), 



or shortly: bp- , F = 0/0, j — c,s. The advantage 
of this transformation is that the term quadratic in the 
bosonic operators in V c ~ c contains now only &0 Cq -type 
operators. The ring Hamiltonian finally reads 



# r ing = £0 ^ E \ n l\ h \jq h F ]q ~ l**M 
Fj 9/0 

- V ?^+,A S p-nf7)W E u._ 1 4-- 



Fba 



£ o E ( ^ + $*SD-(PWrb* ) + ^ooAA 2 



+ 2 E 4|n «' Wq - q { b ® cq + 6 ® c "«) ( & ® c -<? + h ® cq 



9^0 



(18) 



1. Diagonalization 



The Hamiltonian (JT5J) can be diagonalized by a Bogoli- 
ubov transformation. More details can be found in the 
Appendix. As a result we get 



^ing = E E £ F 3 q a F ]q a F M + £ ° E 



A/" 2 

JX Fba 



Fba 



Fj q¥=o 
e 4> J2 sgn(6)AA FbCT + ^ Wx>AA 2 - ^/V. (19) 



Fba 



The definition of the energies ep, as well as the relation 
between the new (ap ■ ) and the old (bp- ) operators can 
be found in the Appendix. The degrees of freedom which 
are affected by the interaction are those related to the 
a ®e<j-operators. Therefore, these are called charged de- 
grees of freedom, whereas the indices (Qcq), (®sq) and 



(Qsq) denote the neutral modes. In total, there are one 
charged and three neutral modes. 

In order to investigate the combined effects of gate volt- 
age and magnetic flux, it is convenient to introduce ap- 
propriate linear combinations of the number counting op- 
erator J^Fra ■ 



N = A/© c :=^AW 

Fba 

J = J®c '= J~] sgn(b)J\f F b<T, 



(20) 



Fba 



being the total particle and total current operators and 
■Afs>s ■= y^sgn(a)J\f Fb(T , J^s := ^2 sgn(ba)Af F ba , 

Fba Fba 

M Qc ■= E sgn(F)A/> b(T , J ec := ^2 sgn(Fb)Af P ba, 

Fba Fba 

Af es := ^sgn(Fcr)A r F 6CT , J Qs := ^ sgn(Fba)Af Fba . 



Fba 



Fba 



Expressing H Ilng in terms of these newly defined opera- 
tors yields 



H,, 



H, 



ring.b 



H Ti 



EE 

Fj 9#0 



s~. at a-. 

Fjq Fjq Fjq 



E (A^1 J +J| J )+ const, 

Fj^®c 



£0 
16 



(21) 



Mg 



and kd 



where W = 1 + MmL , k u =- w 

co ' PS Woo+eo/8 

The first line corresponds to the bosonic part of the 
Hamiltonian, H Y - , and the second and third lines to 

the fermionic part H T - f . An eigenbasis of H Ylng is thus 
formed by the states 

n 4^ #)f= : {#•*)}' (22) 



. Fi,«/o, 



m #i ,! 



where N,0) has no bosonic excitations. In general, any 

number of electrons can be distributed in many different 
ways on the branches {Fba), which is described by the 
set of vectors 

N~K rh N -K rl,N- Ko i h N- Ka ii). (23) 

For each value of N one has to include the states which 
contain the bosonic excitations of the interacting elec- 
trons, where the the parameter mp, counts the number 

of excitations in the channel Fj with momentum g. 



IV. LINEAR TRANSPORT 

We have now all the ingredients to evaluate the linear 
transport characteristics of the interacting SWNT ring 
by use of the Kubo formula. 



A. Conductance formula 

The current operator at lead L is defined as the rate of 

change of particles at the contact L, i.e., 1^ = —eN^ — 
-%{H,N L ] = -f[H TL ,N L ] where -e is the electron 
charge. It yields 



^ = -J E/ d M^(r^t(f)$ CTL (r1 - h.c). (24) 
Analogously the current operator at the right lead is 
f x = f E / d 3 r{T R {r)^ a {r)<$> aR (r) - h.c). (25) 

To proceed, we use relation ([S]), relating the 3D operator 
^(r) to the slowly varying ID one ipFba(x). We assume 
that the latter does not change significantly in the tun- 
neling region. This yields for the current operators 



h 




Fba q 

d 3 rT L (r)<p* Fb (f)^- jL (f) ifj Fba . (x L )c q ~ < 



ah 



ie 

J 



E E T LFbq^ Fba (xL)Cq-aL ~ k.C. 



(26) 



Fba q 



Ir = jJ2J2[ TRFb ^FbA x R) c fcL~h.c.\, (27) 

Fba q \ / 



where xl and x R are in the middle of the respective 
tunneling regions. The above expressions for the cur- 
rent operators can in turn be used to evaluate the DC- 
conductance of the interacting SWNT ring in terms of 
the Kubo formula 



G = lim Re-^ 

n^o nil 



dt'jW-V{[i L {t), I R (t')]) eq , (28) 



with /l(^) and I R (t) being current operators in the in- 
teraction representation, where the interaction is rep- 
resented by the Hamiltonian H ext . Hence, ( ) cq de- 
scribes the average with respect to the equilibrium den- 
sity operator p cq := Z- x e~^ Iin « + ^ T+ ^ L+ ^ ', with /3 
the inverse temperature and Z the partition function. 
It is convenient to introduce the linear susceptibility 



X (t) = -%6(t){[I L (t), /fl(0)]) cq in terms of which J28 
assumes the compact form 



G 



- lim i . 



(29) 



Here x(O) is the Fourier transform of the response func- 
tion x(i). In 4 the linear susceptibility at imaginary times 



x(ri-T 2 )=-(T T [I L (n),I R (T 2 )}) 



cqi 



(30) 



where T T indicates time-ordering, has been evaluated for 
a spinless Luttinger liquid. Generalizing 4 - to the multi- 
channel situation represented by a 3D toroidal SWNT, an 
expression for the linear susceptibility can be obtained to 
lowest non vanishing order in the tunneling matrix ele- 
ments T a p bq : A summary of the calculation, where we 
made use of the explicit form of the nanotube Bloch wave 
functions, can be found in the subsection below. It de- 
livers for the conductance the remarkably simple result 



G 




(31) 



(U!,XL ~Xr) 



where np(w) is the Fermi function and G r pl a (w, Xl — Xr), 
the retarded Green's function for the interacting elec- 
trons on the SWNT ring, is entirely determined by the 
slowly varying ID part of the electron operator ([6]) as it 
is the Fourier transform of 

G$b(t,x-x') = -~0(i)({^ bCT (x,£),4 6 >',O)}). 

(32) 

As the number J\f of electrons in the dot can vary, 
( ) indicates the thermal average with respect to the 
grandcanonical equilibrium density matrix /9 r ing,AA — 

^rTns ex P _/3H ' ing - I n contrast, knowledge of the 3D char- 
acter of the SWNT Bloch wave functions is encapsulated 
in the tunneling functions Q R and <&l, see Eq. (|42p be- 
low. We notice that Eq. (|3"Tj) is not trivial, as it pre- 
dicts the absence of interference between Green's func- 
tions with different indices, in contrast to the case of a 
strictly one-dimensional ring considered in 4 and the for- 
mula given ir>2&. This result has its origin in the very 
strongly localized character of the p z orbitals, and on the 
fact that we assumed extended contacts coupling equally 
to both sublattices of the underlying graphcnc structure. 
For temperatures further than ksT from a resonance (|31[) 
further simplifies to 



G 



l<&d 2 |$7 



l> 2 



E 

Fba 



GTbai^L 



■xr)? 



(33) 



B. Proof of the conductance formula (|31|) 

In this subsection, which the hurried reader can skip, 
the linear susceptibility x(£l) entering the conductance 



formula (129)) is obtained by first evaluating the Fourier function method to find an expression for the imagi- 



transform x(ifl n ) of the imaginary time response func- 
tion x( T ) an d successive analytic continuation: x(fi) = 
linijQ^-tsi x(iQ n ). Specifically, we use the generating 



J 



nary time susceptibility at lowest order in the tunneling 
couplings 4 . It reads 



x(n -t 2 ) = -rj 



h.i 



dr 



hii 



dr ' z2 z2 z2 z2 z2z2 T LF 1 b 1 qTi F2b2ir T RFab3(r T^ Fib ^, x 



Fibi F 2 b 2 F363 F464 era' qq' 



<\{T T ^F lfclCT 0, X L )lpF 2 b 2 cr(ri, X L )lp F3b3(J ,(T' , X^lpFibicr' (T2, Xr) )G La (q,T - Tl)GRa{^ ,t' - r 2 ) + 

- ( T r V'F 1 6 1(T (T,2;L)'0F 2 6 2CT (Ti,a; L )^ 3b3(T ,(r2,X i? )V'F464O-'(T',a;i?.) )G L a{$,T-Tl)G Rcr (( l ','P i -T') + 

- ( T r ^F 1 b 1 A T ^^ X L)^F 2 b 2 a{T,X L )ll) F3b3(jl (T',X R )^ Fibil y l (T2,XR) )G Ler (q,Tl - t)G R(T (^ , / - T 2 ) + 

+ ( T r -0F 1 6 1 <r( r i' :z; i) , /'F 2 b 2CT (T,x L )^ 3b3(T ,(T2,a;/ ? )-0F4fc4<T'( r '^fl) )Gx CT (g, n - t)G Rl7 (<f, r 2 - r') | , 

I 



(34) 



where G aa (q,T — t') = — (T' t (c ( jV tq (t)c-, (t'))) q denotes Though the four particles correlator can in principle 

the Green's function for the free electrons of lead a. In be evaluated using the bosonisation approach, see e.g£, 

(fM)) ( } a describes the average with respect to the equi- following 4 - we factorize it as 
librium density operator Z~ 1 e~ l3Ha of the lead a. 



{Tr^ Flbl A T ^ X L)i>F 2 b 2 a{Tl,X L )%l; F3b3a \T\xR)%k Fibl a'(T2,XR)) 
~ 5F 1 F 2 5b 1 b 2 5F 3 Fjb 3 b i (T T l(j F2 b 2 a(Tl, X L )^ FlbllJ {T, X L )) (lpFib 4 a' (t 2 , Xr)^^, (t* , Xr)) 
- SF 2 Fjb 2 b 3 S F4Fl 6 b2b3 6aa' (T T ^F 2 r 2 cr(Tl, X L )%jJ F3b3a ,(T\ Xr)) (^p 4bia , (t 2 , X R )%/j Fibia (T, X L )) . 



(35) 



In this approximation multiple interference processes are 
neglected. The conductance is thus expressed in terms of 
the single particle Green's functions 

G Fb „{T-T',x a -Xp) := (T T ip Fbr7 (T,x a )tl; Fbry (T' \xp)). 

(36) 
In frequency space this yields 



then follows as 



x(*n») = -|iftfl E *>»."»- 



h 4 h(3 

X /,/,/,/ TLFbqTlp'b' q T RF'b' q'TRFbq' 
Fb F'b' cr qqi 

X G F i b > a (iuJ n , XL - XR)G Fbr7 (iuj n - i^n, Xr - X L ) 

x [G La (q,iuj n -iri n ) ~G La (q,iuj n )} 

x [GRa(q ,iuj' n ) - GR a (q ,iu' n + iVl n )). (37) 



As a further step we perform the summation over the 
fermionic frequencies and carry out the analytical contin- 
uation iil n — > il to find x(£l). From (|29p the conductance 



G 



hh 4 



/ , / , / , / TLFbqTt F 'b'q T RF'b'q'T RFbS , 
Fb F'b' a qqi 

duG^l a (u>, X L - X R )G F % a (u!, XR - X L ) 



dni 



MY 



I — I 4ImG£ (g, w)ImGg (g» , 



duo I 



where the superscript "ret/adv" refer to re- 
tarded/advanced Green's functions, respectively. 
Such expression can be simplified further by replacing 
the summation over the ^-values with an integral over 
energies: Y.q = f depL{e) J2q le , where p L (e) is the 
density of states of lead L. Analogous procedure holds 
for the right lead. Recalling the explicit expression of 



8 



the tunneling function T aFbq ~ one then finds 

C: ™£££ / '*/*(*)£ / '* WOE 

1 .•* .* 



Fb F'6' 



d 3 r 2 Tl(r 2 )ipF'b'{r2) 
d 3 r 3 T R (f 3 )ip* F , b ,(r 3 ) 
d 3 r A T^{r4)ip Fb {r 4 



V L 

1 



1 



(38) 



IVr 



-*«ii 



^G r f t(w,ii - x fl )G^ v (w,x fl - «z,) 



x I -^|M j 47 r 2 ,5(c - s/h)S(u - e'/ft), 

where we made use of the expression for the spec- 
tral density of a free electron gas A aa (q, cu) := 
— 2hxi(G T ^(q, lj)) — 2ttS(uj — e a $/H), and we assumed 
that the leads wave functions are well approximated by 
plane waves. A crucial simplification follows now from 
the observation that the p z orbitals, entering the prod- 
ucts <Ppb(ri)(pFb(r4,)i see (Jl]), are strongly localized at 
the carbon lattice sites. In contrast, the other position- 
dependent functions entering (|38|) are slowly varying on 
the scale of the extension of the localized p z orbitals. 
Hence, we can replace the latter with delta-functions cen- 
tered at the position R + f p of the carbon atoms. The 
summation over the wave numbers q associated with an 
energy e yields 

1 






Jq\e(R+T p -R"-T p „) 



(39) 



4tt sin(Q| e |-R + f p - R" - f v „ |) 



% 



\R- 



R" 



A7tS r,r" 5 p 



where the latter equation means that we assume the 
lead wave vectors q\ eF at the Fermi level to be larger 
than 1/ao, with ao the nearest neighbour distance on the 
graphene lattice. This is satisfied e.g. for conventional 
gold leads. We thus obtain 



G 



( 27 r) 2 e 2 C 2 L 2 vvvv 

~H~¥lV T ~ 2^ 1^2^ 2^ JFbpjFbp'fF'b'pfF'b'p' 
L FF i w 



RR' 

fOC 



,i(F-F')(R*-R' x ) 



PP' 

T l (R + t p 



\T R (R' - 



VJ 



du)p L {hh))p R {huj)\ 



driF(oj) 



xGp ba (u},x L - x R )Gp% /a {uj,XR - x L ), 



(40) 



where the constant C results from the integration over 
the p z orbitals. Since we assume an extended tunnel- 
ing region, the fast oscillating terms with F = —F' are 



supposed to cancel. Furthermore, we assume that both 
sublattices are equally coupled to the contacts, such that 
the sum over R (R 1 ) should give approximately the same 
results for p — 1 and p — 2. Hence, we can separate the 
sum over p and p' from the rest and exploit the relation 



/ , JFbpfFb'p — Sb.b 



Introducing now the transmission function 



R,p 



o)\ 2 Pl(sf), 



(41) 



(42) 



and the analogous one, <J?_r, for the right lead we obtain 
Eq. (|3"Tj) when the lead density of states is approximated 
with its value at the Fermi level (wide band approxima- 
tion) . 



V. THE RETARDED GREEN'S FUNCTION 



As the essential information on the conductance is con- 
tained in G r F l a , the main task lies in obtaining an ex- 
pression for this retarded Green's function in the case of 
a toroidal SWNT. 



A. Time dependence 

First we need to determine the time-dependence of the 
ID electron operator introduced in ©. To this extent, 
it is convenient to make use of the so called bosonization 
identity, where the ID electron operator is written as a 
product of a fermionic operator, the Klein factor, and a 
bosonic field operator. A derivation of this identity can 
be found e.g. in~. It yields 

^Fba(x) = ?7FbCTe i T(A^sg„(6)+0)z e ^ CT (*) e ^ F!jCT (^ 

where the Klein factor rjpba lowers the number of elec- 
trons with given quantum numbers Fba by one, and 



<t>UAx) = iY,HTl e ~ iSMb)9Xb F 



i(b)? C 



-aq/2 



(44) 



q>0 



The infinitesimally small number a > in (|4"4")) serves to 
avoid divergences due to the extension of the linearised 
dispersion relation. Thus the time evolution is deter- 
mined by the evolution of the Klein factors ri Fbcr {t) and 
by that of the bosonic fields 4>Fba{t) (144|) . Let us then 
have a closer look at the Hamiltonian (|2"Tj) . It is made 
up of a part depending on the bosonic operators and one 
depending on the number counting operators. These two 
types of operators commute with each other. Therefore, 
the time evolution of the operators rjFba{t) and 4>Fba{x, t) 
can be treated separately. The bosonic part ff r ing,b yields 
the simple time evolution a- (t) = e~^ ep ii f a~ , and 

Fjq x ' Fjq 1 



analogously for the operator a„ . (£) . This can be related 
to the time evolution of the operators OFo-g appearing in 
the fields <\>Fro- The corresponding calculation for the 
Klein factors is a bit longer and yields 



VFba(t) = 

as well as 



YFb<r(t)r)Fba(0), 



f]Fba 

(45) 



--■ y^WvUM- 



VFba 

(46) 



B. Separation of the correlation function 

Let us now turn to the averaging process. As the num- 
ber of electrons J\f in the dot can vary we have to work 
in the grand canonical ensemble, summing over all values 
of N '. The correlation function 

(46>'>ohWM)) = 

= Ti-[p lineM ^ Fb<T (x',0)^ Fba (x,t)} , (47) 

is a thermal expectation value with respect to the grand- 
canonical density matrix operator p r ing,Af — ^■e~ /3Hring , 

and the trace is over the many body states N,m), cf. 

(f2"2"| . where N and rh define the fcrmionic and bosonic 
configurations, respectively. Accordingly, the trace is 



Tr{...} 



£ £ T(N,m 



,v 






N,m). (48) 



{N}m 



In order to make use of the unitarity of the Klein fac- 
tors, we exploit the fact that the above correlator only 
depends on the time difference. Thus, with the help of 
the definitions ((43J) , J44|) and £46]) we get 



^k>',o)VWM)) - (^ Fba (x',-t)^ Fb(r (x,o) i 

/5ri„ g e~ "^ {x ' -~ t] e-^^ {x ' -*5 



= EEE(^> 



where, with £«ng 1 f(iV) being the eigenvalues of .Hring.f) 
(^ Fbir ^,-t)^ Fba (x,0)) =:£ E ^e-^.^ 

and 



6 ^h^ 



E°" 



e ^ rmg.b 



-#lr(, CT (»',-*) -i<j} Fb<J (x' ,-t) ic/> FbiT (x,0) p i4> Fba (xfl) 



x e ^Fi-. 



0,m 

(51) 



In the bosonic part we only need to take care of the state 
0, rh ) , as the bosonic excitation spectrum is the same 
for all charge states. The partition functions read 

,-/3H xins<t 



z f = E E (#.° 

M {N}u 



EC 



, 7Ti e 



-/3H ri „ g ,b 



0,m\ 



A,0), (52) 
(53) 



In order to evaluate the Green's function we also need 
the correlator (jpFb t r(x,t)'ip Fb(T (x',0)). As above, it can 
be factorized in a fermionic and a bosonic part. It holds, 

ypFba^'^-^Fbaix^)^ = \lp Fbo{x,t)l}) FbiJ (x' ,0yj . 

(54) 

For the fermionic part we find 

VWM)4 b >',0)) =:E E ^-e-'W 
x rF&^e-^C^^M+^K*'-*). (55) 

C. The bosonic part of the Green's function 

The evaluation of the bosonic part of the Green's func- 
tion is lengthy but standard. Following e.gi^i we obtain 



v y} " I t , \ 1 

x e -^^- ssn(6)+ ^V F U-^k CT (0)rWo) \<W* .-*)^M)) 6 = y^- 



(56) 



x e 



¥(^p 



gn(6)- 



3 <^pi«,(a'>0) e »^Fl.<r( a: ' ) 



A^,mV(49) 



1 v-^ e 



-cnq 



4 ^^ ?7 
g>0 q 



As the bosonic and the fermionic operators commute, 
the correlation function can now be separated into a 
fermionic and a bosonic part: 

(4-6>',-*)twm)) = 



X;Q^(-t,sgn(6)(a;'-x)) 
+ 2Cl cqq [Q cq (~t, (x' - x)) + Q cq (-t,-(x' ~ x))}] } . 



where we introduced the correlation function 



Q lq (t,x) :— isinl^2 q X \ 



(57) 



1 — cos( —£ qx 



coth 



ps^ q 
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The parameters C^ entering in the Bogoliubov trans- 
formation which diagonalizes the bosonic part of the 
Hamiltonian, cf. (fT9|) . are defined in the Appendix A. 
To proceed in the analytic evaluation of the bosonic part 
of the Green's function we perform the following two ap- 
proximations. The first one consists in linearising the 
energies s cq which belong to the diagonalised Hamilto- 
nian ((IS]): 



V F 

z cq = hv c \q\ = h-—\q\, 



(58) 



w 



2 I sw <io-<io 



where we introduced the charge velocity v c . 
The coupling constant K c is given by K c — 

.-1/2 

, go = 27T/L, and is determined by 
the interaction for small values of the momentum q. In 
case of a repulsive interaction K c is always smaller than 
1, yielding v c > v F , s cq > s 0q and s c = e /K c > e . 
The second approximation refers to the parameter 
C^cqq- We already assumed in section [Til CI that, due to 
the long range character of the interaction, a momentum 
cut-off q c can be introduced. For high values of the 
momentum q, the parameter C^ cqq , which depends on 
the interaction, has to vanish. We assume an exponential 



decay in g-space, C\ 



(Bcqq 



- C 



®cq a l 



-?/9e. The relation 
between C^ cqo and K c is found with the help of the 
transformation parameters a q and (3 q : 



1 



£()</ 



1 



it i, ' ^"olN,_ 



2s 



2eoq 



For q = q they read: a qo = \^jK c /2, /3 q . 
l/(2-\Z2K c ), leading to the expression 



SWNT Green's function 



1- z[ a+ — ,0 

1c 



E 



E 



i 



e -/«W e - 



■i2g.(N Fh „Bga{b)+4>)(x'-x) 



M {NU 



\f 



OO 

E 

k,l=0 



F k ,i(b,x' - x) 



(61) 



hw - £ u (N Fba 1 N 1 ^ g ,4>)+ irj 

F^(b,x'-x) 

ttoJ - £-k-i{N Fba + 1,N + l,Hg,<j>) + ii] 
where r\ is a positive infinitesimal and the energy 

£-k-i(N Fb(T + l,Af+l,fj, g ,4>) =-M ff + (62) 

e (N Fba + k+- + 0sgn(&)) + W 00 (Af+ -) + sj 

is the energy needed to add a particle in the branch (Fba) 
to a system with M electrons together with a difference 
of k neutral and I charged bosonic excitations. Likewise, 
-£ k i(N Fbcr ,N Cl /j, g , <p) is the energy gained by removing a 
particle on branch (Fba) from a system with Af electrons 
plus bosonic excitations. 

VI. RESULTS FOR THE CONDUCTANCE 

We are now finally at a stage in which we can discuss 
some features of the conductance of a toroidal SWNT, 
as they follow from the conductance formula 1(5T|) to- 
gether with the expression for the Fourier transform of 
the Green's function (l6"T1). 



C 



_ I [ | _ _ ( i^l/2 7^—1/2^ /rn\ A. Summation over the relevant configurations 



®cq g 4 I £ 



With these two approximations we can evaluate the 
bosonic correlation function at low temperatures, ksT <C 
so (in this regime coth(/3eo/2) ~ 1). We find, see Eq. (|B7|) 
in Appendix [Bl 



^PbA^^-^Fbaix^)) = 
I b 



1 — z \ a 



J2J2F k!l (b,x'-x)e^ okt e^ 



ejt 



J $>cq Q 



(60) 



fc=0 1=0 



where z(a, x) — e a e e l z x . 



D. The Fourier transform of the Green's function 



Despite its intricate form, see (|B7|) . the time depen- 
dence of the bosonic part of the Green's function is triv- 
ial. Gathering then together bosonic and fermionic con- 
tributions we obtain for the Fourier transform of the 



The expression (J61J) is quite intricate, as it requires a 
summation over all the fermionic and bosonic configura- 
tions. A noticeable simplification, however, occurs in the 
low temperature regime s$ > kgT, which we shall ad- 
dress in the following. In this regime in fact we need to 
consider ground state configurations only of the SWNT 
ring, i.e., we can neglect bosonic and fermionic excita- 
tions. Neglecting bosonic excitations means to set k = 
I = in the sum in ((61]) , such that the only effect of the 
bosonic part of the Green's function is in the interaction- 

r 2 



dependent prefactor 



1 — z ( a 



-.0 



With a = 



a + q c 1 it yields a conductance suppression proportional 

to (d/L) 2C ® c " = (a/L) (K - +K - 1 - 2 ^ 4 with respect to the 
non-interacting case. To find the relevant fermionic con- 
figurations we observe that at low temperatures the con- 
tribution to the trace comes from those values of the total 
electron number J\f and of the total current J which min- 
imize the energy £ , r ing,f(A r ). This energy is provided by 
the fermionic part of the Hamiltonian (pH]) and is deter- 
mined by the applied gate voltage and magnetic field. 
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EkflA 



FIG. 3: Ground states resonances for a ring with two 
fermionic species b = R,L: Along the lines the energies are 
degenerate and transport is there maximal. Continuous lines 
correspond to the interaction parameter W = 2, dotted lines 
to the non-interacting case, W = 1. 



1. Spinless case 

To clarify the ideas, let us first consider a ring with 
only two species of spinless electrons characterized by 
b = R,L. Then, H r ing,t simplifies to 



H, 



ring,f 



£o 
4 



W(N - k^) 2 + (J - k^f 



(63) 



where now k u 



Pa 



-20 and FT = 1 



2W 00 

*»9 - Woo+£o/2> ""f ~ " Y """ " ~ i_r eo 

hold. Figure [3] shows which ground state (N,J) is occu- 
pied in the (Au, fe„ )-plane. The continuous lines are ob- 
tained by setting W = 2, whereas the dotted ones belong 
to the non- interacting case W = 1. Along such lines the 
energies of neighboring ground states are degenerate. At 
zero temperature these are the conductance resonances. 
Figure [4] shows the ground state energies E{k ll ,7V, ^7) 
for k$ = 0. For values of k^ around zero the ring is in 
the state (AT = 0, J = 0). Raising the external gate volt- 
age above k^ = 0.75 leads to a transition to a charged 
state with one electron. Above k^ = 1.25 the state 
(J\f = 2, J = 0) is the one with the lowest energy. We 
observe that along k^ = for the states with an odd 
number of electrons there is also a degeneracy between 
the states with the same charge but different current di- 
rections. 




FIG. 4: Ground state energies E(k llg ,jV,J') for k,/, = and 
interaction parameter W = 2. Notice that the parabolas cor- 
responding to an odd electron number are doubly degenerate. 



2. SWNT ring 



Let us now turn back to the SWNT ring. Due to the 
presence of the spin and orbital degrees of freedom, the 
fermionic part of the Hamiltonian is quadratic in all of 
the variables Np, and Jpp cf. (f2i"j) . However, it re- 
mains gapped only with respect to the total charge and 
current operators. Let us then introduce for given /j, g and 
4> the integer numbers A^ c and J c which minimize .Ering.f ■ 
Outside resonance, only one couple (N c , J c ) contributes 
to the grandcanonical sum in (|6ip . Nearby resonance 
we have to include states which differ by one electron 
and with corresponding current configurations differing 
by plus/minus one. Moreover, in the low temperature 
regime /c^T < £o the quadratic forms in H rmg [ impose 
that we only have to take into account those configura- 
tions {N}, where the number of electrons in each con- 
figuration (F'b'a 1 ) differs at most by one (i.e., we ne- 
glect " fermionic excitations" ) . From now on we measure 
the number of excess electrons or holes from the charge 
neutrality point where N = 0. Thus, N F'b'a' can only 
assume the values m and m + 1, where m > above 
and m < below the charge neutrality point. This 
means that N Cl obtained by summing over the occu- 
pation of all the branches (F'b'a 1 ), can be written as 
N c = 8m + l, m G Z, I = 0, 1, ... 7. Notice that this 
implies an eight-electron periodicity as a function of the 
gate voltage. In contrast, as J c counts the occupation 
difference between the orbital branches, b — R, L, it can 
only take values between —4 and +4. We can thus sim- 
plify the trace operation considerably by summing only 
over the dominant configurations in (|6ip . These are the 
configurations where N c — > 2V C + 1 with Npba — Tn, and 
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N c + 1 -> iV c with 7V_f &ct = m + 1. We obtain 



£$>,*'-*) = 



1 - : ( cH ,0 



e -^ ri „ E ,f(^l«c + i + e -/5 S r i„ g ,f(^liVc + i) e -i^ S gn(&)(x'-x) 



e "' 9E r 1 n g ,f( JV |iVc) + e "/3S ing ,f(JV|jV c + l) 

e -i^(msgn(h) + 0)(a:'-2;) 



fiu — £oo(m + 1,N C + l,/i g ,0) + irj 



(64) 



Recalling that ^ ring , f (^V|tv c+ i ) - £ ri n g ,f(^k+i) = 
£oo (w + 1, iV c + 1, /i 9 , ^>) = £q, we obtain 



1- z[ a-\ ,0 

1c 



\G%l(u,x'-x)\ 2 = 

cos 2 (Ax) + sin 2 (Ax) tanh 2 (^fa) 
\hu - £ 00 (m + 1,N C + l,/j, g ,<f>) +ir}\' 2 



zC.< 



(65) 



where Ax := ^(msgn(6)(x' 



When Eq. 



is 

inserted in Eq. (J3TJ) the conductance diverges, as a con- 
sequence of the presence of the infinitesimal convergence 
factor irj rather than of a finite line-width^S. Neverthe- 
less, Eq. (|6"5"j) allows already to discuss the main quali- 
tative features of the conductance, namely the resonance 
pattern as a function of the applied gate voltage and 
magnetic field. 

To be definite, conductance resonances occur at a de- 
generacy for the removal of either a counter-clockwise or 
a clockwise propagating electron. Moreover, as in the 
spinless case, double resonances occur for special values 
of the applied magnetic flux and gate voltage, such that 
clockwise and anti-clockwise propagating modes are si- 
multaneously at resonance. As the ring energy is inde- 
pendent of the spin and Fermi number degrees of free- 
dom, a double resonance is uniquely fixed by N c , J c . Let 
us consider e.g. N c — 0, J c — and N c + 1 = 1, J c ± 1. 
Then it only exists a configuration with N c — 0, while 
there are four equivalent configurations with N c = +1, 
J c = 1 and four with N c — 1, J c = — 1 contributing to 
the double resonance. 

The resonant pattern is easily obtained by approximat- 
ing the derivative of the Fermi function in Eq. (|3ip with 
a delta function^ such that G ex ^ Fba \G r pl a (0, x' — x)| 2 , 
and then evaluating the poles of (f6"5)) as a function of k^ 
and k^. Upon recalling that N c = 8m + I, with m € Z 
and I = 0, 1, ..7, we find 



fc Mg = 8m 



4 
W 



m 



III! 



so 



8Woo, k^ 

—~fl ~ 7^ s S n ( & ) • ( 66 ) 

e Q W W 



Notice that this implies a distance between adjacent res- 
onant lines, with I and I + 1, of Ak u = ^a. The 
corresponding results are shown in figure [5] We choose a 
value of W — 2 for the plot. The solutions for the non- 
interacting system are shown with dotted lines. In the 
non-interacting case the ring is occupied by a multiple of 




FIG. 5: Conductance resonances of the carbon nanotube ring 
as a function of k lJig and k$. The interaction is chosen to be 
W = 2. The lines with different slopes correspond to an excess 
of electrons propagating in opposite directions (counterclock- 
wise and clockwise). Dashed lines are for the non interacting 
case. 



four electrons. From one square to the next other four 
electrons can be accomodated at one time, as there is no 
charging energy. In the presence of interactions, however, 
the picture is entirely different. The former squares are 
deformed to hexagons which only touch at their vertical 
boundary. There, the two ground states are degenerate. 
The conductance is zero here. In the parallel stripes in 
figure [5] we indicated one possible fermionic configura- 
tion only. There is no difference in the energy if we add 
the first electron to a spin-down state instead of a spin- 
up one as we did. The choice of F is arbitrary as well. 
The main difference to the non-interacting case is that, 
similar to the spinless case, for certain values of fc M no 
conductance resonances can be achieved. 

Finally, Eq. (|cj5|) predicts a periodic pattern of the 
conductance as a function of the distance between left 
and right leads, with maxima when the distance between 
leads is xr — xl = L/2, cf. Fig. Q] 



VII. CONCLUSIONS 

In this work the linear conductance of an interacting 
toroidal SWNT weakly coupled to leads and pierced by 
a magnetic field has been analyzed. We focussed on low 
temperatures ksT < e$, where £o is the mean level spac- 
ing. In this regime Coulomb blockade effects are known 
to crucially affect the conductance of a straight SWNT 
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as a function of an applied gate voltag o 15 i 16 . Likewise 
Coulomb interactions strongly influence the conductance 
characteristics of the ring. Due to charging effects, the 
electrons can enter the ring only one by one and not 
up to four at the time as in the non-interacting case. 
While in straight tubes the conductance pattern has a 
four-electron periodicity, for toroidal SWNTs we predict 
a periodicity of eight electrons. The conductance exhibits 
Coulomb resonances, and the resonance pattern is a func- 
tion of the applied magnetic field <j> and gate voltage /i ff . 
In contrast to the non-interacting case, where the param- 
eters <f> and /ig can always be tuned to match a resonance 
condition, in the interacting system this is no longer pos- 
sible and we find that at certain values of [i g a window 
opens where no conductance resonances can be obtained. 
Interference effects are manifested in the occurrence of 
double resonances where electron trajectories propagat- 
ing clockwise and anticlockwise are degenerate. Finally, 
the conductance as a function of the distance between 
left and right leads is maximal when the distance be- 
tween leads is half of the circumference length, i.e., when 
clockwise and anti-clockwise propagating electrons have 
to cover the same source to drain distance. 

The conductance of toroidal SWTNs has already 
been measured in the experiments^ 3 ., though not in the 
Coulomb blockade regime. Thus, we believe that the pos- 
sible verification of our predictions is within the reach of 
present experiments. 



APPENDIX A: BOGOLIUBOV 
TRANSFORMATION 

The non-diagonal part of (TT5]) is given by 

H nd = eoJ2\ n i\ b i>c q b®cq + ^^2^\n q \W q ,- 



q^O 



q^O 



x [b^cqb^c-q + 6© C q6 ffiC9 + 6 ffic _ g & ffic _ g + 6 ffic _ g 6 ffiC9 



kt 



I-,t 



hi 



hi 



(Al) 

The goal of the Bogoliubov transformation is the deriva- 
tion of new operators, in terms of which iJ n d is diagonal. 
In general, we can write them as a linear combination of 
the old bosonic operators b^ cq and b^ cq : 

a lcq = Y.( U 1l' b icq>-Vqq'b ( Bc q '), (A2) 

a ®cq = ^(-Vqq'b^q, +U qq ,b (Sc q>). (A3) 



Of course, these operators have to fulfill the bosonic com- 
mutation relations: 



a ec9> a e C( j<] : 


— Oqq' , 




(A4) 


Qffiaj, a® C( j'] : 


= [ a (Bcq, a (Bcq'\ 


= 0. 


(A5) 



Additionally, in order to ensure the diagonal form of H nd 
in terms of the new operators, the conditions 



[^c,a® 



,t 



[H c , aq 



-e cq a mq (A6) 



have to be fulfilled. Note that also new energies e cq have 
been introduced. The calculation of the transformation 
matrices U and V follows closely the guidelines in 33 . The 
result is 



U 




V = 



where eo q = £o\n q \ and the transformed energies e cq are 
given by 



eo g yl + 8W g ,_g/eo. 



Using the abbreviations a q = t> 



gQcj 



and f3 q = \ 



(A9) 



'2>-\i H 



the new operators read, assuming q > 0: 

a <S>cq = \ a q ' Pq){"®c—q ~>~ "®ci}i 

+ (fiq-a q )(b @c - q + b @cq ), 

a®cq = {Pq-OLq){b^ c _ q +b^ cq ) 

+ (a g +/3 g )(6® c _ ? -(-6® cg ), 



a 1 

u ©c- 

a ec- 
The inverse of 1 

°©cg 


q 

i 

q = 
:he 

+ 

+ 

+ 
_i_ 


= {a q +l3 q ){b^ c _ q -b\ cq ) 

- (ot q - (i q )(b 9c - q -b mq ), 

-- (a?-/3 9 )(&© c _ ? -&© cg ) 

- {a q + P q ){b( S c~ q -b( Bcq ). 

transformation is 

-{<*q + Pq)(a®c-q- a ®cq) 


8a q (3 q 
(a q -/3g)(a©c- ? + a© cg ) 


b($cq 


8a q P q 

(a q — Pg)(a 0c _ 9 + a fficg ) 


8a q [3 q 
-(a q + /3 q )(a ec _ ? - a© C g) 


°®c-q 


8a q f3 q 
("g+/3g)(o©c-g + °©cg) 


8a q (3 q 
-(a q - /3 g )(a ffic _ g - a® cg ) 


b($c-q 


8a q (3 q 

— {Uq — Pq){a^ c _q — d^cg) 


8a q f3 q 
(a q +/3 g )(a ec _ g + a ecg ) 



8a Q l3, 



qHq 
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The inverse transformation can also be written in a more 
concise way: 



h Fjq ~ Z2 \ S F 3qq ' a Fjq' + C Fjqq' a 



t 
Fjq 



where we have for the spin related coefficients S and C: 

b($sqq' = t>Qsqq' = &Qcqq' = Oqq' > 
^®sqq' = ^Qsqq' = ^Qcqq' = 0, 

while for the charge related ones, we have to take care of 
the sign of q. Let us choose q > 0: 

o _ r aq+ Pq _ , «g+/?g 

„ _ i- Qq " Mg , r a q ~ Pq 



Be— gg' — "gg' 



c f 



Sao' 



a q + 0q , a q + Pq 
8a q f3 q q q 8a q j3 q 

0^ q p q Oi q p q 

" Ua—a' ~~ 



c - qq ~ " qq 8a q p q uq - q 8a q p q 



The final result for the diagonalised Hamiltonian is then 

Ml 



Pj q^O Fra 

+ £ ^sgn(r)AA Frff + il¥ooAA 2 , (A10) 

Fra 

where we introduced the energies 

£ ®sq — SQsq = £ Qcq = £0g an d £($cq = Ecq- (All) 

APPENDIX B: THE BOSONIC CORRELATOR 

We evaluate in this appendix explicitly the bosonic 
part of the Green's function, Eq. (|56|) . in the so called 
Luttingcr limit, where i) we assume a linear disper- 
sion: San = £oq/K c , with coupling constant K c — 



cq 



\J\ + 8W qo - qo /eo; ii) we assume an exponential decay 
of the parameter C^ cqq entering the Bogoliubov trans- 
formation: C^ cqq = C^ )Cqa e~ q / q ", with q c a momentum 
cut-off. Let us then start from Eq. ([56]) : 



(<,>',-*)VW^o)) = - — ^ 

\ / 6 l — e a 



(Bl) 



x exp 



g>0 



Z^ QFjq 
.Fj 



(t,sgn(b)(x' - x)) 



+ 2Cl coo [Q cq (t, (x' - x)) + Q cq (t,-(x' - ar))]] } ; 



dcqqV°icq\^, \^ ■" } } T ^cq\ 

where we introduced the correlation function 



*°Cyq{ti X) '.- 



i sm 



£ry q T 



1 — COS 



qx 

(s^ q t 



(B2) 



qx 



coth 



m 



In the low temperature regime we can approximate 
coth(/3eog/2) and coth(/3e cq /2) by 1. We perform now 



the Luttinger approximation and recall that q = 2nn/ L. 
Rewriting the sin and cos functions as exponentials and 
remembering the formula 



pun 

Y / — = -\n(l-e a ), \e a \<l, 



(B3) 



n>0 

we arrive at 



-"ec9 



(^ Fb a(x',-t)iJ Fba (x,0)) 

- r* 

l-z(a-\ ,0^ 

— zla, vpt + sgn(b)(x' — x) 
1 — z(a,v c t + sga(b)(x' — x)j 

1 — z ( a H , v c t + (x 1 — x) ) 

\ q c /. 

1 — z I a H , v c t — (x' — x) ) 

_ ^ c . 

where we have defined the function 



-3/4 
-1/4 



2.r' 2 



_ 1 ,^2 

2°ec? 



£(a,x) 



— p(-a+tai)3f 



(B4) 



(B5) 



Because the functions entering (|B4[) are periodic in time, 
we can write the bosonic correlation function as a product 
of Fourier series. To this extent we use the relation 



(i-*r-E 



fc=0 



r(Ar-7) fe i i 1 

*Mr' |z|<1 ' 



which yields the final result 



(B6) 



CD'-'QO 



\*pFbA x ' ',-t)lpFba{x,0)) = 1-Z I <X-\ ,0 

x ^r fe la,--,sgn(b)(x'-x) \e^ £okt 

fc=o V ' / 

00 

x ^s;(6,x' -x)e^ £jt . (B7) 



i=0 

Here we have defined 



''A: 



( a 7 x ) ■= r ( fc -7) fc ^(- a+ „) 

1 ' 7 ' j - fcir(- 7 ) 



s/(&,o;'-.t) := ^ r ni I a, --,sgn(6)(x' - a;) J 

j^ecgo' ( X ~ X ) J 



I— Til 

x ^ r» 2 I m 

n 2 =0 
/— n± — ri2 

x > r„ 3 I o 

«3=0 



n 3 =0 

where T is the gamma function 



1 1 

1 1 



' 2 ® c ?o' ' X *) ) ' 
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